Filtering of the absolute value of photon-number difference 
for two-mode macroscopic quantum states 



O 

<N 

bjo: 
< 



43 

a*. 



> 

O 

00 

o 



1 



M. Stobiriska, 1 ' 2 F. Toppel, 1 ' 3 P. Sekatski, 4 A. Buraczewski, 5 

M. Zukowski, 6 ' 7 M. V. Chekhova, 1 G. Leuchs, 1 ' 3 and N. Gisin 4 

1 Max Planck Institute for the Science of Light, Erlangen, Germany 
2 Institute for Theoretical Physics II, Erlangen-Niirnberg University, Erlangen, Germany 
3 Institute for Optics, Information and Photonics, 
Erlangen-Niirnberg University, Erlangen, Germany 
* Group of Applied Physics, University of Geneva, Geneva, Switzerland 
5 Faculty of Electronics and Information Technology, 
Warsaw University of Technology, Warsaw, Poland 
6 Institute for Theoretical Physics and Astrophysics, University of Gdansk, Gdansk, Poland 
7 Hefei National Laboratory for Physical Sciences at Microscale and Department of Modern Physics, 
University of Science and Technology of China, Hefei, Anhui, China 

We propose a filter that selects two-mode high number Fock states whose photon-number differ- 
ence exceeds a certain value. Such a filter is important for the engineering of macroscopic quantum 
states of light and for the control of bright light beams. It improves distinguishability of some 
states but preserves macroscopic superpositions. We suggest an operational implementation of such 
a filter. To this end we employ an interference effect in which two input ports of a beamsplitter 
are entered by highly populated multi-photon states. If the photon numbers at the two input ports 
are almost equal, the photons are highly asymmetrically distributed between the two outputs. If 
photon numbers differ a lot, the exit ports get nearly equal photon occupation numbers. 



The set of quantum states that are efficiently produced 
is limited. Nevertheless, with quantum state engineering 
one can modify or enhance certain properties of exper- 
imentally accessible states. In particular, measurement 
induced state operations allow one to filter out states 
of required properties. In the case of intensity mea- 
surement, of considerable help are threshold detectors, 
selecting Fock states or superpositions with sufficiently 
high population numbers. Examples of low- threshold de- 
tectors are known: single photon on-off detectors and 
human eyes [l], Q . Such detectors can be utilized in feed- 
forward techniques based on tap measurements, in which 
effectively POVM measurements [3[ are performed lead- 
ing to quantum operations. As a result, it is possible to 
block light of unwanted properties (too low or too high 
intensity). More complicated filters for Fock states are 
based on the interference effects 0, 0] • 

A more challenging task is to construct a filter that 
would select states with certain properties without de- 
stroying quantum superpositions. This is especially im- 
portant for macroscopic superpositions of Schrodinger- 
cat type, which suffer from low distinguishability and 
are easily destroyed by inefficient detection Such 
macroscopic superpositions combine quantum properties 
with macroscopic population and could enable efficient 
light-matter coupling, and therefore are of much interest 
for various protocols in quantum information technology, 
such as quantum memory jjj [To| , quantum key distribu- 
tion [ll| , macroscopic Bell inequality tests [12|, ll3| , and 



ter (MDF). It filters out two- mode states of light whose 
mode populations differ by more than a certain threshold. 
Its optical implementation will be based on the interfer- 
ence of multi-photon states entering a beamsplitter, each 
mode via a different input port. Theoretically, the aimed 
filter should perform the projection operation 
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quantum metrology [14] 



Here, we propose a blueprint of such a device which 
is further called a modulus of intensity difference fil- 



where |fc, Z) is a two- mode Fock state. Further on, with- 
out the loss of generality, we will assume that the two 
modes of the initial state are distinguished by orthogonal 
polarizations, while all other characteristics are identical. 
For 5 t h > the filter acts as "quantum scissors" [l5j. It 
cuts out those terms of an expression of the state in the 
Fock basis for which the occupation difference is below 
the threshold and preserves those where the difference is 
above it. 

We would like to comment on the two key features 
of such a filter. It estimates the absolute value of the 
difference instead of the difference itself. This has an ad- 
vantage. Since the spectrum of the operator Vs th is fully 
degenerate, the filter does not provide any information on 
which mode was more populated. Further on we discuss 
the action of a MDF on a macroscopic qubit encoded in 
polarization (see formula fl3J). We show that it is not able 
to distinguish the states and filters them fairly. Secondly, 
the estimation is performed in a "yes" -"no" manner: the 
exact value of the modulus is never measured and thus, 
this information does not exist. This property is of great 
importance for all quantum protocols that require state 
engineering preserving the superposition. This is why we 
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FIG. 1: (Color online) Comparison of two filtering techniques: 
absolute difference (MDF) (a) and orthogonality filter (b). 



call this device a filter. 

These two features are the main difference between the 
MDF and the orthogonality filter (OF) based on intensity 
difference measurements [16] . OF is the basic element in 
the setups performing measurement induced operations 
on the macroscopic polarization states [12| . It splits the 
state in two orthogonal polarizations on a PBS and mea- 
sures their intensities I and I± by photon multipliers 
(PM) separately. This is a von Neumann measurement 
TTki = l)(k, l\ and it projects the state on a single Fock 
state. Next, the electronic signals from PMs are passed to 
the OF and electronic "dichotomization" is applied. The 
intensity difference is calculated by taking the signals dif- 
ference and it is compared with a threshold value. The 
+1 (-1) value is signaled if / — 1± > 5th (I± — I > 5th), 
and otherwise. Thus, this is a comparative threshold 
measurement. Contrary to MDF, such a scheme destroys 
superpositions. Measurement induced operations on the 
multiphoton states based on the OF device allow only 
for efficient state discrimination, not filtering. Thus, OF 
is not suitable for preselection strategies in Bell tests, of 
the kind suggested in [l2|. There, polarization analy- 
sis on the reflected beam leads to state identification in 
the transmitted beam and breaks the entanglement. The 
action of MDF and OF is compared in Fig. HJ The for- 
mer projects onto Si and S2 areas simultaneously and 
all superposition terms belonging to them are preserved. 
The latter combined with photon multipliers projects the 
state on a Fock state either in S\ or S2 (red or blue dot 
in the figure). 

Let us analyze the action of the MDF on specific 
"macroscopic qubits," produced by the optimal phase co- 
variant quantum doners in the process of phase sensitive 



coherent parametric amplification [H, [l6|, [TtJ , 
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where (p and (p 1 - are the two orthogonal polar- 
izations from the equatorial plane of the Poincare 
sphere, the real- valued probability amplitude 7^ = 
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FIG. 2: Distribution p$ computed for g — 1.87 and filtering 
threshold 5 th = 200 with 0% (a), 50% (b), 90% (c) losses. 



cosh<T 2 ((tanh0)/2) i+j ^(1 + 2i)\(2j)\/i\/ j\, g is the 
parametric gain. Please also note that, due to the dif- 
ferent parity of occupation in the Fock basis, |$) and 
|3>j_) are orthogonal and thus form a basis for a macro- 
scopic qubit. In a recent experiment [16|, each pulse con- 
tained up to 4sinh(g) ~ 10 4 photons on the average. 
However, in such a photon number regime, detectors are 
not single photon resolving, but distinguish counts vary- 
ing by ±150 photons in the best case [17]. This makes 
the macroqubits hardly distinguishable through direct 
detection [16|. Distinguishability may be quantified in 
terms of the photon distribution p<$>(k,l) = | (k, l\<f>) | 2 
and p® ± (k,l) giving the probabilities of finding simul- 
taneously k photons in cp and I in cp 1 - 
00 

P*(k,l)= ^2 7^Mi+i^,2j, p<s> ± (k,l)=p<$>(l,k), (3) 

i,j=0 

where the Kronecker delta 5 a ^ equals 1 if a = b or 



3 



1 

0.8 
0.6 
0.4 
0.2 






= 


Sib. 


= 40 


— ■ — 5th 


= 200 


5th 


= 400 



0.2 0.4 0.6 0.8 1 
R 




FIG. 3: Distinguishability v evaluated for the gain g — 1.87 

and different threshold values 5 th as a function of the losses FIG - 4: An implementation of modulus of intensity difference 
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filter. The description is in the main text. 



otherwise. Distinguishability is given by 

v = Pi Sl) - = Pi Sl) - Pt\ (4) 

where P$ S ^ = ^2 k leS , p&(k,l) is the probability of find- 
ing |$) in Si of Fig. [Hand P^ Sl) +P^ S2) = 1. It increases 
if |$) (|3>_l)) starts to occupy mostly one of Si regions, 
e.g. Si (S2), with increasing 5 t h- Fully distinguishable 
(indistinguishable) states have v = 1 (v = 0). 

Our analysis is performed taking into account possi- 
ble losses, which are modeled by a BS with reflectivity 
R giving the loss. The p$ distributions evaluated for 
g = 1.87 and 5 t h — 200 in the presence of 0%, 50% and 
90% of losses are depicted in Fig. [2j The filtering results 
in cutting out a stripe \pl5th wide located symmetrically 
along the k = I line. The state |3>) occupies two disjoint 
regions of space: the top (Si) and bottom (S2) triangles, 
but increasing the threshold reduces the contribution of 
p<$> in 62: the peak value goes down originally from 10 ~ 2 
to 2 • 10" 4 for 5 t h = 200 and R = 0. At the same time the 
distribution peak in Si increases from 10 -2 to 4 • 10 -2 . 
Similar behavior is observed for higher gains. Behavior of 
p<$> ± is identical but mirror reflected. Thus, distinguisha- 
bility increases. The particle loss results in shifting the 
distribution towards the origin of the coordinates, i.e. 
vacuum state. The distribution peaks become smooth 
and symmetric. The edges along the threshold lines are 
blurred and the bigger losses are, the smaller width of the 
gap is. It disappears completely for 90% of losses. With 
increasing losses the height of the upper and left peak 
first, drops due to the fact that the distribution "melts" 
and next, increases because the total probability over the 
shrinking area has to be 1. 

For states (j2j) we computed their distinguishability v 
for several filtering thresholds 5 t h as a function of losses, 
see Fig. 02 If no filtering is applied, the distinguishability 
v = 0.64 is independent of the gain value and losses, but 
drops quickly to if R > 0.9. If 5 t h increases, v increases 



as well and approaches unity with reasonable probability, 
e.g. v = 0.96 with the probability of success p s = 10 -4 . 
Obviously, for R = 1 we get v = 0. 

A physical implementation of an MDF can be realized 
by a tap measurement, interference of Fock states on a 
beamsplitter and the feed forward technique, see Fig. HI 
Let us consider a superposition of two mode polariza- 
tion Fock states. The tap measurement is realized by a 
highly biased BS with low reflectivity r, e.g. r = 10%, 
which reflects only a small fraction of the state leav- 
ing the transmitted part almost unchanged. The anal- 
ysis of probability distribution for a BS shows that, for 
highly populated input with super-Poissonian statistics, 
such BS should not reflect more than 20% of the input 
intensity to preserve most of the character of the ini- 
tial state in the transmitted beam. Under such condi- 
tions, an estimation of the population difference for the 
reflected beam will give us an estimation for the trans- 
mitted beam. One should employ here the "worst" case 
assumption. The reflected part will be filtered via an 
interference and measurement process, and the measure- 
ment result will determine a feed forward action on the 
transmitted beam, which will be opening or closing of a 
shutter. 

Let us describe in detail the filtering setup. The anni- 
hilation operators of the two polarization modes of the re- 
flected beam will be denoted aS Qjf , CLrp-L . The correspond- 
ing intensities are given by X r = a£a r , X r ± = a [ r± a r ±. 
They enter the input ports of the polarizing beamsplit- 
ter (PBS). We shall assume that the PBS is oriented in 
linear polarization HV basis with the output operators 
being an = 1/ V^(a r + )•> a v — 1/ V2(a r ± — a r ) . The 
two exit beams, polarized H and V, propagate to a pair of 
detectors, which measure the intensities. If at a given mo- 
ment the respective detectors' readings give photon occu- 
pation numbers K and L, this is equivalent to a collapse 

of the state into \^) out = \K,L) = l/\/K\L\ a ] R K a ] y L \0) . 
This state, according to the aforementioned mode trans- 
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FIG. 5: Distribution of the population difference p(S) in the 
input state with s = 200 photons conditioned on the mea- 
surement of K = 0, L = 200 (a), K = 60, L = 140 (b), 
K — L — 100 (c) at the PBS output. 



formation at our PBS can be rewritten in the terms of 
the input modes as the following superposition [l8[ 
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V(Pl+P2)W( K + L -Pl-P2V- 
\pi + P2, K + L - pi - P2)r,r±- ■ 



(5) 



The analysis of the superposition shows that if the out- 
put photon numbers K, L are measured and the cases 
where they differ little are selected, then, with a high 
probability, the input photon numbers would differ by 
more than a certain threshold. Vice versa, if K and L 
differ a lot this implies small differences in the occupation 
numbers of the input modes, see Fig. [5j where character- 
istic examples are given. In this way we can filter out the 
states for which the photon-number difference exceeds a 
certain threshold value. Note that the states in terms of 
input modes occupation numbers are not Fock ones, but 
superpositions of them. 

Fig. [5] shows the probability distribution p(S) of the 
population difference in the input state S = K + L — 
2pi — 2p 2 for given K and L and the total photon num- 



ber s — 200. We also calculated the probabilities that 
| S | was below or above a given 5 t h (vertical dashed lines 
5th — 30). The filter works probabilistically. For any 
outcome K and L all values of S are possible, but not 
equally probable. For the majority of the cases, for a 
reasonable threshold the probability of high population 
difference is higher than the probability of the low one. 
If the probability is satisfactory e.g. above 90%, as is for 
K = L = 100, the feed forward signal opens a shutter in 
the transmitted beam of the initial BS, Fig.HJ Otherwise, 
we reject it and the shutter remains closed. 

The proposed scheme of modulus of intensity difference 
filter is feasible. This is a threshold measurement which 
projects the input state on a certain multi-dimensional 
subspace of the Hilbert space, thus enabling quantum su- 
perpositions and further state processing. The filter em- 
ploys the interference of multi-photon states on a beam- 
splitter and it works for any two-mode polarization state 
with super-Poissonian statistics. This property is impor- 
tant due to the properties of the tapping BS. Realization 
of such measurements with commonly used photon detec- 
tors is quite challenging, however properties of the MDF 
are worth the possible effort. The filter could be applied 
in the engineering of macroscopic states of light. In the 
case of macro-qubits it circumvents the problem of in- 
efficient detection, and improves their distinguishability, 
thus making them good quantum information carriers for 
quantum information protocols and quantum metrology. 
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